New bounds on lepton flavor violating decays of vector mesons and the Z° boson 



Thomas Gutsche 1 , Juan C. Helo 2 , Sergey Kovalenko 2 , Valery E. Lyubovitskij 1 " 

1 Institut fur Theoretische Physik, Universitat Tubingen, 

Kepler Center for Astro and Particle Physics, 
Auf der Morgenstelle 14, D-72076 Tubingen, Germany 



o 

(N 
C 



Oh- 

> 

m 



o 



X 



2 Departamento de Fisica, 
Universidad Tecnica Federico Santa Maria, 
Centro Cientifico-Tecnologico de Valparaiso 
Casilla 110-V, Valparaiso, Valparaiso, Chile 

(Dated: January 14, 2013) 

We give an estimate for the upper bounds on rates of lepton flavor violating (LFV) decays M — > 
fj^e* of vector mesons M = p°,ui,<p,J/ijj,T and the Z° boson in a model independent way, analyzing 
the corresponding lowest dimension effective operators. These operators also contribute to nuclear 
H — e-conversion. Based on this observation and using the existing experimental limits on this 
LFV nuclear process, we show that the studied two-body LFV decays of vector bosons are strongly 
suppressed independent on the explicit realization of new physics. The upper limits on the rates of 
some of these decays are significantly more stringent than similar limits known in the literature. In 
view of these results experimental observation of the two-body LFV decays of vector bosons looks 
presently unrealistic. 
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It is well known that lepton flavor violation (LFV) is 
strongly suppressed in the Standard Model (SM) by a 
very small neutrino mass. Therefore, the observation of 
any LFV process would be a signal of physics beyond 
the SM. Various avenues can be devised to study these 
issue. When elaborating on a strategy for searches of 
LFV these processes should be considered, which have 
the best prospects for discovery. This includes both 
the prospects for experimental identification of the LFV 
events and possible theoretical limitations on the corre- 
sponding rates. Latter considerations may, for example, 
deal with model independent relations between different 
processes, some of which are already strongly limited by 
experimental data. 

In the present paper we study from this point of view 
LFV decays of vector mesons and the Z°-boson 



at the VEPP-2M e + e collider. They give a model in- 
dependent upper limit on the — > e/i branching fraction 
of 



Br( 



e/i) < 2 x 10" 



(2) 



Also, there exist experimental limits for the e/i decay 
mode of J/ip and of the Z° boson 

Br(J/ip -> e/i) < 1.1 x 10~ 6 , 
Br(Z° -> e/i) < 1.7 x 10" 6 

and /it decay mode of T [2[ : 

Br(T -> fir) < 6.0 x 10~ 6 . 



(3) 
(4) 



M 



with M = p°,U),4>, J/ip,r,Z°. (1) 



The abundant production of vector mesons and Z°- 
bosons in current experiments naturally suggests to 
search for their two-body decays in the /t ± e T final state, 
which is rather convenient for event identification. Re- 
cently the SND Collaboration at the BINP (Novosi- 
birk) [l[ reported on the search for the LFV process 
e+e~ — > e/i in the energy region y/s — 984 — 1060 MeV 



In the near future this list may be extended by the results 
of other experimental collaborations. However, a natural 
question, which arises in this context, touches upon the 
prospects of this category of searches in view of possible 
theoretical limitations on the rates of these LFV decays. 

In the literature there already exist stringent limits of 
this sort. For example, unitarity relations between the 
vector boson LFV decays given in ([1]) and the pure lep- 
tonic LFV decay fi — > 3e have been exploited in Ref. 0] ■ 
From the existing experimental bounds on the latter pro- 
cess the following stringent bounds were deduced Q 
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Br(0-> e/i) < 4 x 10~ 17 , 
Br( J/i> -> e/i) < 4 x 10~ 13 . 
Br(T 4e/i)<2x 10~ 9 , 
Br(Z° -> e/i) < 5 x 10~ 13 . 



(5) 
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In the present article we approach the LFV decays ([T]) 
from another point of view relating these processes to 
nuclear p — e conversion which is tightly constrained ex- 
perimentally [J[ . 

In Refs. [a, Q we studied the LFV process of nuclear 
p" — e~ conversion in the framework of an effective La- 
grangian approach without referring to any specific real- 
ization of physics beyond the SM responsible for LFV. We 
examined the impact of specific hadronization prescrip- 
tions on new physics contributions to nuclear p~ — e~ 
conversion and stressed the importance of vector and 
scalar meson exchange between lepton and nucleon cur- 
rents. In particular, we derived limits on various LFV 
couplings of vector mesons to the p — e current using ex- 
isting experimental data on p~ — e~ conversion in nuclei. 
The purpose of the present paper is to use these limits to 
set upper bounds on the rates of the vector boson LFV 
decays given in ([1]). In Ref. [7j we already indicated in 
this framework upper limits on the rates of the LFV de- 
cays of p°, u and 4> mesons. Here we extend our analysis 
to J ftp, X and the Z° boson and compare our results with 
exiting experimental bounds ©, ([3]) and the theoretical 
predictions of Ref. Q . 

The contribution of vector bosons to p~ — e~ conver- 
sion in nuclei is shown in Fig. 1. In this diagram the 
upper vertex corresponds to the LFV interactions of vec- 
tor bosons M — p°,u>, <j>, J/i/>, T, Z° with e,fi given by 
the following model-independent Lagrangian [UJij]: 



MM 

'"nff 



M"(e^ +h. c .) . a,', 



Here the £y A are effective vector and axial couplings of 
a vector boson M to the LFV lepton currents — ^-liil 1 
and j£ = ejfj,j5fi. The possible effect of additional non- 
minimal derivative couplings of vector bosons to the LFV 
lepton current will be considered later. 

The lower vertex of the diagram in Fig. 1 is described 
by the nucleon- vector boson Lagrangian Q: 



1 



Cmnn = -Nj^N y g 



M 



(7) 



where g MNN are effective couplings. In this Lagrangian 
we neglected the derivative terms which are irrelevant for 
coherent p~ — e~ conversion. The Lag rangian jOmnn is 
an extension of the conventional nucleon-vector meson 
Lagrangian (il-fToj. In the case of the light p°, uj and cj> 
mesons we use values for g MjjL ^ which are taken from an 
updated dispersive analysis 0, 0j| 



4-0 , 9b _ 



41.8, g, NN =-0.24. 



(8) 



In addition we need an estimate for these effective cou- 
plings involving J/ip, T and Z°. The couplings g J/ilNN 
and g rNN can be extracted from data [2| on J/ip — > NN 
and T — > NN decays: 



T(J/ip -> pp) = 202 ± 9 eV , 
r(T ^pp) < 27 ± 1 cV. 



(9) 



JV N 

FIG. 1: Contribution of intermediate vector bosons to nuclear 
/i _ — e~ conversion. 

The two-body decay rate for M — > pp, where M = J/i/j 
or T, is given in terms of the effective coupling constants 

9mnn 

by: 



r(M -> pp) 



12tt 



■ m M 1 



4™^ 



m 



A I 



2m 



N 



(10) 



Using the central value of T(J/ip — > pp) and the upper 
limit for T(T — > pp) Q we get 



w 1.6 x 10~ 3 , 
< 3.3 x 10" 4 . 



(11) 
(12) 



For our purposes the upper bound, as far the latter case, 
coupling is not sufficient. As it will be seen later from 
Eq. (|20|) . for the derivation of upper bounds on the ef- 
fective couplings £,v.A entering in ^ we need a definite 
estimate for the value of g rNN - This can be done on the 
basis of a QCD analysis of exclusive processes of heavy 
quarkonia as M — > pp [l2l [iH ] . The corresponding tran- 
sition amplitudes are generated by three-gluon annihi- 
lation between the heavy and the light quarks. In this 
approach the following expression for the coupling gMNN 
of a heavy quarkonium state M of mass uim with the nu- 
cleon was derived as 



fh 



ga s {m M 



2 s /M 



(13) 



' M 



Here g = 95.4 GeV 4 represents the loop integral over 
nucleon wave functions. Using data for the leptonic decay 
constants = 416.4 MeV and / T = 715.5 MeV @ we 
get: 



9 J/iPNN 
9rNN 



0.14 al{m 2 JM ), 
0.9 x 10~ 3 a 3 (m|) 



(14) 
(15) 



coincides with 



Above value (fT4l for the coupling g J/qflIS , 
the value (jTTJ) extracted from the experimental data 
for the strong coupling constant a s (m^) = 0.226, 
which is in a good agreement with the world average of 
a a (m 2 ji^) ~ 0.26 [lij]. This exercise can also be regarded 



3 



as a consistency check of the approaches evaluated in 
Refs. [H|,[ll|]. Similarly, using the central value of the lat- 



15] extracted from 



est result for a s {m 2 c ) = 0.184 
radiative T decays we find from Eq. (| 15[) the value of 



+0.015 
-0.014 



9i 



= 5.6 x 10" 



(16) 



This result is significantly smaller than the upper bo- 
und (1121) set by experiment (for a discussion see also 
Ref. (3). 

The coupling of a Z° boson to nucleons is well known 
(see, for instance, [13] )■ Since the axial nucleon current 
does not contribute to the dominant coherent channel of 
nuclear fi — e conversion we only need the coupling to the 
vector nucleon current, which in view of Eq. ([7]) we denote 
as gzNN- Neglecting a possible but small contribution of 
str ang e and heavy sea quarks in the nucleon it is given 

as \m 



9, 



9 



2 cos 9w 



(1 -4 sin 2 6 W ) « 0.31 



(17) 



Here we used the following values of the SM parame- 
ters: M w = 80.399 GeV, sin 2 9 W = 0.2322, G F = 
1.16637 x 10~ 5 GcV~ 2 . 

Starting from the Lagrangian £^ = C 1 ^ + Cmnn 
of Eqs. ([6]) and (JT)) it is straightforward to derive the 
contribution of the diagram in Fig. 1 to the total fi~ — e~ 
conversion branching ratio [f| . To the leading order in the 
non-relativistic reduction the coherent fi~ — e~ conversion 
branching ratio takes the form [18| 



rycoh Q 



PeE e 



— , (18) 
2tt — > capture) 

where p e , E e are 3-momentum and energy of the outgoing 
electron (for details see Ref. [EEl); r(/z _ — > capture) is 
the total rate of the ordinary muon capture reaction. The 
factor Q in Eq. (18]) has the form dlHI 

Q = \{M P +M n )a% + (M P -M n )a^ v \ 2 + 
+ \{M V + M n )a { % + {M P -M n )af v \ 2 . (19) 

It contains the nuclear matrix elements A4 Pi n which have 
been calculated numerically in Refs. [l9j for various nu- 
clei. Here we consider pT — e~ conversion in 48 Ti studied 
by the SINDRUM II Collaboration 0. For this nucleus 
we have M p « 0.104, M n « 0.127. The Q factor also 
contains the LFV lepton- nucleon parameters avv,AV- 
For the contribution of the vector boson-exchange dia- 
gram in Fig. 1 these coefficients are expressed in terms 
of the LFV couplings ^/ of Eq. © as [5-7] 



In Ref. we extracted upper limits on the couplings 

a *y from the experimental upper bounds on fi~ — e~ 
conversion in 48 Ti reported by the SINDRUM II Collab- 
oration |4J. These limits can be translated into bounds 
on the LFV couplings £yj{ . Assuming that no acci- 
dental cancellations occur between the different terms 
in (|20|) and using the values of the vector boson-nucleon 
couplings from Eq. JSJ, (fTTj). (fTH|> and <[TTJ) we get for 
non-derivative couplings 

^<3.6xl0" 12 , < 3.6 x 10~ 14 , 

it < 1.0 x 10~ n , g+ < 1-4 x 10~ 8 , (21) 

£j<3.5xl0~ 5 , £f ° < 6.4 x 10~ 8 . 

Note that the Lagrangian ([6]) also governs the LFV decay 
M — > efi of vector bosons. Thus, using the limits from 
Eq. (|2"Tj) we can set upper bounds on the rates of these 
two-body decays. Their branching ratios are given by: 



Br(M ->• efi) 



127rr 



M 



m M 1 



M ) i 



(22) 



where Tm = TUp/mu and is the total decay width of 
boson M. Here we neglect the electron mass. With the 
limits set by Eq. (|21[) we get the following upper limits 
on the branching ratios of the vector boson LFV decays 



(p° -> efi) < 3.5 x 10~ 24 , 
(w -> efi) < 6.2 x 10~ 27 , 
{(p^efi) < 1.3 x 10~ 21 , 
(J/tp^efi) <3.5x 10~ 13 . 
(T -> efi) < 3.9 x 10~ 6 , 
Br(Z° -> efi) < 8.0 x 10" 



Br 
Br 
Br 
Br 
Br 



(23) 



-15 



The limit for the J ftp — > efi mode is compatible with the 
corresponding number extracted from /i — > 3e as done in 
Ref. [3[ and shown in Eq. ([5]). For the cases of the <f> — >• e/i 
and Z° efi decay modes we obtain significantly more 
stringent upper limits than in Ref. Q while for T — > efi 
our limit considerably weaker. The existing experimental 
upper bounds on some of these rates listed in @ and Q 
are by far much weaker than the limits set by theory, 
both as in Ref. [3j and as discussed here. 

Finally, let us comment on the effect of non-minimal 
derivative couplings of vector bosons to the LFV lepton 
current in the upper vertex of the diagram in Fig. 1. 
One could imagine a situation where the minimal non- 
derivative couplings in the effective Lagrangian JSJ) are 
substituted by the following derivative couplings [3j, [5[ : 



* (3) 



2 m 2 + m M 



2 

9 HNN 



(20) 



2 ^ m 

H 



where a = V, A; H — w, <j>, J/tp, T, Z°; thm and m„ the 
vector boson and muon masses, respectively. 



MM 



t¥H„ +^jl)M^+h.o., (24) 



where M^ v = d^M v - d v M^ is the stress tensor. Here 
the are effective tensor and pseudotensor couplings 

of the bosons M to the LFV lepton currents jj„ = 

ea^ufi and = ecr^„ 75/i. As was noted in Ref. 
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the derivative couplings (|24|) would lead to significant 
weakening of the bounds in ([5]). This happens be- 
cause the contribution of the virtual vector bosons to 
fi — > 3e is reduced in comparison to the case of the 
non-derivative couplings ([6]). This is also true for nu- 
clear /i — e-conversion and the derived bounds of (|23[) . 
These bounds would have to be divided by the follow- 
ing factors @: q 2 /m 2 M m 2 l /(2m 2 M ) — 10~ 2 [p°,u;]; 5 x 
10~ 3 [^];5.7 x 10~ 4 [J/#6.1 x 10" 5 [T]. However, even 
with this weakening the limits ()23[) still exclude experi- 
mental observation of the LFV decays of vector mesons 
p°,uj,cf),J/i(},T — > /i ± e T in the near future. The situ- 
ation with derivative coupling of the Z° boson to the 
leptons is different from the case of the mesons and 
was studied in Ref. [2(| ■ This coupling is induced by a 
SUl(2) x LV(l)-invariant effective operator of dimension 
higher than four. Due to the electroweak gauge invari- 
ance couplings both of the Z°-boson and of the photon 
to the LFV lepton current are induced by the same op- 
erator. Therefore, Z° — > c\x can be constrained from 
the existing experimental data on /i — > ey and the elec- 
tron electric dipole moment. In this way a very strin- 
gent bound Br(Z° -> efj,) < 10~ 23 — 1CT 22 was derived 
in Ref. [2(J. This is significantly more stringent than 
both our bound in Eq. (23) and the bound (5;) derived in 
Ref. 0. 

In conclusion, we extracted from the experimental 
bounds on nuclear \i~ — e~ conversion new upper lim- 
its on the LFV couplings of the vector mesons and the 
Z° boson to the e — /i lepton current. Then we applied 
these limits to deduce upper bounds on the branching 
ratios of LFV decays of vector mesons and the Z-boson. 
The obtained upper bounds are shown in Eq. (|23[) . Our 
bounds for the decays p°,uj — > ep are new. The bounds 



for <fi — ¥ e/i are significantly more stringent than the cor- 
responding bounds existing in the literature. This con- 
clusion indicates, in particular, that the nuclear ji — e- 
conversion is more sensitive probe of LFV than the de- 
cays of Z° boson and vector mesons. In the latter case 
this is true at least in the sector of the lepton interac- 
tions with the mesons made of quarks of the first and 
second generation. On the other hand searches for the 
LFV decays of Z° boson and mesons remain an impor- 
tant experimental effort since their observation at the 
rates above the limits ([2^]) would be a manifestation of 
new LFV physics, which does not fit into the present 
analysis. In particular, it may imply a non-trivial mech- 
anism of self cancellation of different terms in (fl"§)) . which 
we considered as unnatural. Then /i — e-conversion rate 
could remain below the experimental bound [4J while 
allowing large rates of the LFV decays of Z° boson and 
vector mesons. 
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